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I. INTRODUCTION 



The newly created Bose-Einstein condensates (BEC) of weakly interacting alkali-metal 
atoms [1] stimulated a large number of theoretical investigations (see recent reviews [2]). 
Most of these works are based on the assumption that the properties of BEC are well de- 
scribed by the Gross-Pitaevskii (GP) mean-field theory [3]. 

Recently, it has become possible to tune atomic scattering length to essentially any 
value, by exploiting Feshbach resonances (FR) [4,5]. A fundamental open problem is how to 
describe the physics of dilute BEC near FR (dilute in the sense that the range of interatomic 
potential, R, is small compared with inter-particle spacing r ) in the regime of a large 
scattering length, a, which we take to be positive. The GP approach fails in the regime of a 
large gas parameter, n \ a | 3 , where n is the particle density. 

The dilute BEC for a large gas parameter regime in one-channel approximation has been 
considered previously in Ref.[6] (for the corresponding problem for a Fermi gas see Ref.[7]). 

In this paper we consider the ground state properties of the dilute homogeneous Bose 
gas near the FR using a multichannel zero-range potential (ZRP) model of FR. 

In section II we describe the lowest order constrained variational (LOCV) method [6,8] 
for the one-channel N-body problem. The calculations for model interaction potentials used 
by Ref. [9] are presented. The description of the multi-channel ZRP model is given in section 
III. Section IV develops the LOCV method for the dilute Bose gas near FR. We conclude 
the paper in section V with a brief summary. 

II. LOWEST ORDER CONSTRAINED VARIATIONAL METHOD 

In a dilute many-body problem in the large gas parameter regime correlations between 
particles are very important. 

The LOCV method [6,8] for the homogeneous N-body system is to assume a Jastrow 
many-body wave function, of the form 



*(fi, r 2 , ...f N ) = f(fi - fj), 



(1) 



i<j 




(2) 



while at large distances / must approach a constant. 
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In the LOCV method [6,8] the boundary conditions for f(r) are 



f(d) = 1, /'(d) = 0, (3) 

the expectation value of the energy is given by 

E/N = 2im\ I* f 2 (r)r 2 dr, (4) 
Jo 



and d is defined by the normalization 



Ann 



[ f(r)r 2 dr = l. (5) 
Jo 



In Ref.[6] for the dilute case (R <C r , where the inter-particle spacing r = (3/(47m)) 1 / 3 ) 
inter-atomic interaction was replaced by the zero-range potential (ZRP) model [9] 

(rf)', 1 



rf lr=0 a 

In this case 



(6) 



dsin(fcr+^ 

/lr<dj rsin(A;d + o)' UJ 



with A;coto~ = —1/a, kdcot(kd + 8) — 1, where A; = \Jm\/h 2 . 

In the small a/r limit, 5 = — fca, d = r and the LOCV result for _E/iV [6] is given by 

— = 2tt n 8 

and is same as that first found by Lenz [10]. 

The ground state energy per particle, E/N, in the low-density regime, na 3 1, can be 



calculated using an expansion in power of y/na 3 

E = ^ an[1 + 1?8 (^3)1/2 + g( ^ _ V3)na 3 [ln(na 3 ) + C] + ...]. (9) 
iv m 15a/7t 3 

The coefficient of the (na 3 ) 3 / 2 term (the second term) was first calculated by Lee, Huang, 
and Yang [11], while the coefficient of the last term was first obtained by Wu [12]. The 
constant C after the logarithm was considered in Ref. [13]. 

The expansion (9) is asymptotic, and it was shown in Ref. [14] that the Lee-Huang- Yang 
(LHY) correction (second term in Eq.(9)) represents a significant improvement on the mean 
field prediction (the first term in Eq.(9)) . 
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In Refs. [15-19] the Lenz-Lee-Huang-Yang (LLHY) expansion (first two terms in expansion 
(9)) has been used to study effects beyond the mean field approximation. In Ref.[18], it was 
found that the correction to the GP results may be as large as 30% in the ground state 
properties of the condensate, when the conditions of the JILA experiment for 85 Rb are 
considered [5]. 

In the large a/r limit, 5 = ir/2, kdtan(kd) = — 1, and the energy is given by [6] 

E n 2//3 Ti 2 

— = 13.33/i 2 = 10.2597 5 , (10) 

N m 2mrQ 

which is very close to the Legett's unpublished variational result 

E 3vr 2 h 2 n 2 / ? > 



11) 



N 2 m 
(this was quoted by Baym [20]). 

To study the validity of the ZRP model we consider an example of the square-well (SW) 
potential with a/r — > oo. The calculated energies per particle, E/N, are presented in Table 
I. The ratio R/tq is typically of the order of 10~ 2 . From Table I, we can see that even for 
°/ r o ~~ *■ °° the ZRP model is a very good approximation (the difference for E/N between 
the ZPR and the SW is less than 1% if R/r « 10~ 2 ). In this case the LOCV results for 
E/N have universal properties that depend on the interatomic potential only through the 
single low-energy parameter a, even for large gas-parameter regime. 

We note that for an attractive potential the atomic BEC is metastable and energy per 
particle can be written as E/N — iT/2. Therefore the real part of energy in LOCV method 
would be reliable if E/N > T/2 [6]. 

In Ref . [14] , the authors have used a diffusion Monte-Carlo method to calculate the lowest- 
energy state of uniform gas of bosons interacting through four different model potentials that 
have all the same scattering length a. 

(i) The hard-sphere potential 

oo (r < a) 



(12) 
(r > a), 



where the diameter of the hard-sphere is equal to the scattering length, 
(ii) Two soft-sphere (SS) potentials 

' V (r < Ro) 



V (ss \r 



(13) 
(r > Ro), 



with R = 5a, V = 0.031543^7 (W) and R = 10a, V = 0.003408ft 2 / (ma 2 ). 
(iii) The hard-core square-well (HCSW) potential 

+oo (r < R c ) 

v {HCSw ){r) \ _ Vo {Rc<r<Ro) (14) 
(r> Rq). 

The parameters of the HCSW potential are R c = a/50, R = a/10 and V = 412. 815a 2 /(ma 2 ). 
The HCSW potential has a two-body bound state with energy -1.13249a 2 /(ma 2 ) [21]. 

Comparison the LOCV results for potentials (12-14) with the available diffusion Monte- 
Carlo (DMC) calculations (Table II, Fig.l and Fig. 2) shows that the LOCV energies in the 
case of R/r 10~ 2 are in very good agreement with the DMC results. 

III. ZERO-RANGE POTENTIAL MODEL OF FESHBACH RESONANCE 

We start from the coupling channel equation 



h 2 

— V 2 + V p (r))i) p (r) + g(r)il> Q (r) = E^ p (r), 



(15) 



h 2 



(-—V 2 + V Q {r) + S)^ Q {r) + g\r)^ p {r) = Eijfi (r ) , 



where S is the energy shift of the closed channel Q with respect to the collision continuum. 
Since potentials V p (r), V®(r) and g(r) have a range typically of the order of interatomic 
potential range or less, we can replace Eq.(15) by 



(--V 2 + V^if) + V l2 ^ Q {r) = E^j p (r), 
m 



(--V 2 + V 22 + £)^ Q {r) + V 21 4> p (r) = E^(r), 
m 



(16) 



where 

47ra 2 8 
V*lf) = -—M^S{f) w r (17) 

is a multichannel generalization of the Huang pseudo-potential [22]. 
Equations (16) can be rewritten as free equations 

f) 2 

VW) = Ei/) p (r), 

m 



(18) 



h 2 

( V 2 + SU Q (r) = Eij Q (r), 

m 



=0 = (M u r^ P + M 12 r^ Q )\ r=0 , 

(19) 

=0 = (M 21 r^ p + M 22 r^)\ r=0 , 



with boundary conditions 

d{r%l) P ) 
dr 

dr 

which is a multi-channel ZRP model [9,23]. 

It is easy to show that the Hamiltonian of the particle moving in the field of multi-channel 
ZPR is Hermitian, that is, the condition 

J [(tfVvVf + (0?rV 2 0^ 3 r = / [(V 2 0f M + (V 2 0?)*0f]rf 3 r (20) 

is valid for any 0f , 0*p and <f>%, <p 2 which satisfy the boundary conditions (19) with energy 
independent constant Hermitian matrix M. 

Indeed, we can transform the volume integrals, Eq.(20), into a surface integrals 

1 = £ / WiY(^W 2 - cf>i(nVMr\dS, (21) 

where i = P,Q, the surface S consists of small sphere around r = and infinite sphere, and 
the normal n is directed outwards from S. The surface integral over infinite sphere equals 
zero, and Eq.(21) becomes 

Using Eq.(19) we can rewrite Eq.(22) as 

i,k 

Thus for the hermitian matrix M (M = M + ) the surface integral / vanishes and hence the 
hermiticity condition, Eq.(20), is satisfied. 

We shall use the following parameterization for the matrix M 

M u = - — , M 12 = M 21 = - A M 22 = -I-, (21) 

where a^g is the background value of the scattering length. The model (21) with 7=1 was 
considered in [23]. 
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Solutions of Eqs.(18,19) can be written as 

il) P (r) oc sin(/cr + 5), 

(22) 

ij Q (f) oc e-v 7 ^, 

with £ = mS/h 2 , k 2 = mE/h 2 . 

Since the energy shift £ can be converted into an external magnetic field B by £ oc B, 
the scattering length a depends on the external magnetic field by the dispersive law [24] 

« = ^(1 + ^^), (23) 

where 

A(5) = (1- y^Xv^B + ^B^ (24) 

a 



4£ ^_ / A 
(25 + A)2' 7 " VA + 2iV ^ 

and A = A(_B ) characterizes the resonance width. 

Two parameters of our model Eq.(21), a^g and /3/7, are completely determined from the 
external magnetic field dependence of the scattering length. For example, in the sodium case 
[25] the width A of the Feshbach resonance (B = 907G) is 1G, therefore a = 1.1013210" 3 1/G, 
/3/7 = 2.3472610" 2 , and a bg = 53au. 

As for parameter 7, it can be determined from the energy dependence of the scattering 
phase shift, 5 

k cot 5 = — (-1 + A =). (28) 

ab 9 A + 2B - ^JAB B - (ka bg (A + 25o)/7) 2 

We note here that mean-field approximations with contact potentials were considered in 
many papers [26-30]. 

IV. DILUTE BOSE GAS NEAR FESHBACH RESONANCE 

For a generalization of the LOCV method for a dilute uniform gas of bosons interacting 
through the two-channel ZRP, Eqs.(18, 19,23), we assume 

V p (r u r 2 , ...f N ) = I] f P (fi - fj), (29) 

i<j 
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where \I> p (fi, r 2 , ...r N ) is the Jastrow wave function, index P denotes the projection onto the 
Hilbert subspace of the incident (atomic) channel, and f p at short distance is solution of 
the Schrodinger equation 



with boundary conditions 



— |W P (r) = \rf p (r), (30) 



\r=0 — r J , 



dr 



l eff 

(31) 



where 



df p 

f(d) = l, |u = o, 

a e ff — a bg (l-\ -, ), (32) 

2B - ^/4B B-(Ka bg (A + 2B )/ 1 y 

and k = VmX/h. The real part of ground-state energy per particle is given by 

E/N = 27mA [ d (f p (r)) 2 r 2 dr, (33) 
Jo 

and d is defined by the normalization 

Aim (\f p (r)fr 2 dr = 1. (34) 
Jo 

We note that the effective scattering length, a e ff is a many-body parameter (it depends on 
E/N), and, for B = B , a e ff does not tend to infinity. 

The calculated energies per particle, E/N, of Na(F = l,mj? = 1) atoms at resonance 
magnetic field (B = B = 907G) for r = 10 2 a bg are compared with the one-channel ap- 
proximation, Eq.(10), in Fig. 3. These comparisons show that for finite values of 7 there are 
significant differences between our results and the approximation of Ref . [6] . Our results are 
much smaller than the Legett's variational estimate [20]. 

To consider, so called, nonuniversal effects [21], i.e. the sensitivity to the parameters of 
the interatomic interactions other than the scattering length, we calculate E/N as a function 
of 7. Table III shows a strong 7 dependence of E/N near the FR. 

However near the FR the atomic BEC is metastable, and the LOCV E/N would be 
reliable if E/N > T/2 [6]. We have extracted the values of T/2 from experimental data of 
Ref. [25]. Using these values of T/2 we have calculated the ratio (T/2) /(E/N) (see Table IV). 
From Table IV we can see that the LOCV results for the real part energy, E/N, is valid for 
the experimental conditions of Ref. [25], since (T/2)/ (E/N) w 10~ 2 . 



Now suppose that a, Eq.(25), near the FR is negative (a bg for the Na(F = l,m F = 1) 
atoms is positive). In one-channel case the uniform condensate for negative a is always me- 
chanically unstable. But the two-channel consideration can lead to stable uniform solution, 
since the many-body parameter a e ff can be positive. Table V illustrates this stabilization 
effect. Although the three-body recombination processes [23,25,29-31] can make it diffi- 
cult to observe this effect experimentally, we note that Ref.[32] considered the possibility 
of suppressing three-body recombinations in a trap. There is a similar case [33] of uniform 
ID gas of N bosons on a ring for which inelastic decay processes, such as the three-body 
recombination, are suppressed in the strongly interacting and intermediate regimes. 
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V. SUMMARY AND CONCLUSION 



In summary we have considered the LOCV method for a dilute uniform gas of bosons. 

Comparison of the LOCV results for potentials Eqs. (12-14) with the available diffusion 
Monte-Carlo (DMC) calculations shows that the LOCV energies in the case of R/r 1CT 2 
are in very good agreement with the DMC results. 

We have generalized the LOCV method [6,8] for dilute uniform gas of bosons near the 
FR, using the multi-channel ZRP model. 

As an example of application, we have considered Na(F = l,m,F = 1) atoms near the 
B = 907G FR. 

At high density na 3 ^> 1, there are significant differences between our results and one- 
channel ZRP [6] for the real part of energy per particle. 

For the case of negative scattering length, we point out the possibility of stabilization of 
the uniform condensate. 
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TABLE I. Energy per particle, E/N , in units of k 2 / (277W 2 ,) vs R/r for the square-well at- 
tractive potential with a = oo. The ZRP model for this case gives E/N = 10.2597a 2 /(2mro) 
[6]. 



R/r 


E/N 


10" 1 


11.1906 


5 x 10~ 2 


10.7143 


io- 2 


10.3502 


5 x KT 3 


10.3057 


IO" 3 


10.2703 



TABLE II. Energy per particle, E/N, in units of 27rh 2 na/m, as a function of R/r for 
soft-sphere potentials. R is the range of potential and r = (3/ (Ann)) 1 / 3 is the atomic 
separation. 



R/r 


E/N 


E/N[14] 


R 


0.081 


1.01961 


1.00427 


5a 


0.161 


1.01951 


1.00427 


10a 


0.174 


1.04292 


1.01382 


5a 


0.347 


1.04288 


1.01302 


10a 


0.374 


1.09599 


1.04167 


5a 


0.748 


1.0973 


1.03689 


10a 


0.806 


1.23315 


1.11011 


5a 
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TABLE III. Energy per particle of Na(F = l,m F = 1) atoms, E/N, multiplied by 
10 6 2mal /h 2 as a function of 7 and B near the B = 907G Feshbach resonance. 



7 


B = 9(JoG 


b = yuo.oCj 


B = 9(Jo.8G 


£5 = 9(J0.9Lr 


O nn/ 1 no/"i 


O nno r (~i 

B = 9(Jo.5Cj 


0.1 


5.147 


6.008 


f* r 

6.835 


(710c 

7.185 


t /i no 

7.498 


n on 

0.889 


0.5 


6.078 


8.939 


15.59 


21.51 


on ro 

29.58 


n nno 

0.998 


1 n 

1.0 


0.128 


n 00 r 

9.225 


18.06 


on r'n 

29.69 


r r on 

55.30 


1 nno 

1.002 


9 n 






1 q nn 


34 Q4 


1 nn 3 


1 nn'? 


3.0 


6.144 


9.319 


19.20 


36.41 


137.9 


1.003 


4.0 


6.144 


9.319 


19.27 


36.41 


168.2 


1.003 


5.0 


6.145 


9.327 


19.30 


37.27 


197.5 


1.003 


6.0 


6.145 


9.328 


19.32 


37.42 


215.5 


1.003 


7.0 


6.145 


9.329 


19.33 


37.52 


231.4 


1.003 


8.0 


6.145 


9.330 


19.34 


37.58 


243.0 


1.003 


9.0 


6.145 


9.330 


19.34 


37.62 


254.5 


1.003 



12 



Table IV. Ratio (T/2)/(E/N) of Na(F = l,m F = 1) atoms vs 7 and B near the B = 
90 7G Feshbach resonance. 



7 


B = 906G 


B = 906.5G 


B = 908 


.5G 


0.1 


2.1 x 10~ 3 


1.8 x 10~ 2 


1.2 x 10" 


-3 


0.5 


1.8 x 10~ 3 


1.2 x 10~ 2 


1.1 x 10" 


-3 


9.0 


1.8 x 10~ 3 


1.2 x 10-2 


1.1 x 10" 


-3 



TABLE V. Energy per particle of Na(F = l,mp = 1) atoms, E/N, multiplied by 
10 6 2mal g /h 2 , effective scattering length, a e ff in units of a^, as a function of 7 at B = 907.1G 
with a(B) = —9ab g . 



7 


E/N 


a eff 


0.1 


8.02 


2.59 


0.5 


49.1 


13.8 


1.0 


146 


31.8 


1.5 


299 


49.5 


2.0 


502 


72.6 
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LU 



2.5 - 



2 - 



1.5 




0.5 - 



10" 7 10" 6 10" 5 0.0001 0.001 0.01 0.1 



a 3 n 



FIG.l. Ground-state energy per particles, E/N,m units of 2%h 2 na/m, vs a 3 n for hard- 
sphere interactions. The LOCV results are shown by full line. The diffusion Monte-Carlo 
calculations [14] are shown as dashed line. 
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LU 



2.5 



1.5 



0.5 - 



10 7 10" 6 10" 5 0.0001 0.001 0.01 0.1 1 



a 3 n 



FIG. 2. Ground-state energy per particles, E/N,in units of 2irh 2 na/m, vs a 3 n for hard- 
core square- well potential (HCSW). Triangles and circles correspond to the LOCV and the 
diffusion Monte-Carlo [14] results, respectively. Lines are drawn to guide the eyes. 



15 



LU 



12 



10 



8 - 



4 - 



2 - 



10 



12 



7 



FIG. 3. Ground-state energy per particles, E/N, of Na(F = l,rriF = + 1) atoms at 907 
G in units of 10 4 2ma^ 9 /^, 2 , as a function of parameter 7, (solid line). Dashed line represent 
approximation of Ref.[6], Eq.(10). (ro = 10 2 ab g ). 
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